BAR4R-CONF-02/08 
SLAC-PUB-9231 
hep-ex/OOOO 



Evidence for the b — > u transition B — ► Dtir and a search for 



The BAEL4R Collaboration 



February 5, 2008 



Abstract 

We report evidence for the b —>■ u transition B — > D~^ir~~ and the results of the search for 
B° — ► D* + 7r~ from a sample of 61.6 x 10 6 T(4f?) decays into -B meson pairs collected with the 
BABAR detector at the PEP II asymmetric e + e~ collider. The observed B° — > D~^ir~ yield has 
a probability of 4.4 x 10~ 4 to be a fluctuation of the background (3.5 a) and we measure the 
branching fraction B(B° -> D+vr") = (3.1 ± 1.0 (stat.)il.O (syst.)) xl0~ 5 . We also set a limit 
B(B° -> £>* + vr-) < 4.3 x 10" 5 at 90% C.L. 
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The measurement of the CP-violating phase of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
|H is an important part of the present scientific program in particle physics. CP violation manifests 
itself as a non-zero area of the unitarity triangle [Q]: while it is sufficient to measure one of the 
angles to demonstrate the existence of CP violation, more than one angle must be measured to 
demonstrate that the CKM mechanism is the correct explanation of this phenomenon. There are 
already several theoretically clean measurements of the angle j3 |§ but there is no such measurement 
of the other two angles (a and 7). A theoretically clean measurement of sin(2/?+7) can be obtained 
from the study of the time evolution of the B° —* D(*^ + ir~ and B° — > D^~ir + ^ decays. This 
measurement requires a knowledge of the ratio between the decay amplitudes of these two processes. 
Unfortunately B° — > D^ + tt~ decays, which are Cabibbo suppressed, cannot be efficiently isolated 
from B° — > D^ + tt~ decays. Since the amplitude of the decay B° — > L>(*) + 7r~ can be related 
to the decay B° — ► Z?i*^ + 7r _ [|j, the measurement of sin(2/3 + 7) will require the knowledge of 
£>(B° — ► Z?i*^ + 7r - ). Decays of this type have also been proposed to be used for the measurement 
of |^ufe/V^b| ||]. In this paper we present the results of the measurement of the branching fractions 
of the decays B° — ► Di*^ + ir^ . 

The data were collected in the years 2000-2001 with the BABAR detector at the PEP-II asym- 
metric e + (3.1GeV) - e~(9GeV) storage ring ||. Since the BABAR detector is described in detail 
elsewhere @, only a brief description is given here. Surrounding the beam-pipe is a Silicon Ver- 
tex Tracker (SVT), which provides precise measurements of the trajectories of charged particles 
as they leave the e + e~ interaction point. Outside of the SVT, a 40-layer drift chamber (DCH) 
allows measurements of track momenta in a 1.5 T magnetic field as well as energy- loss measure- 
ments, which contribute to charged particle identification. Surrounding the DCH is a Detector of 
Internally Reflected Cherenkov radiation (DIRC), which provides charged hadron identification. 
Outside of the DIRC is a CsI(Tl) electromagnetic calorimeter (EMC) that is used to detect pho- 
tons, provide electron identification and reconstruct neutral hadrons. The EMC is surrounded by 
a superconducting coil, which creates the magnetic field for momentum measurements. Outside of 
the coil, the flux return is instrumented with resistive plate chambers interspersed with iron (IFR) 
for the identification of muons and long-lived neutral hadrons. We use the GEANT H software to 
simulate interactions of particles traversing the BABAR detector, taking into account the varying 
detector conditions and beam backgrounds. 

We select hadronic events with a minimum of three reconstructed charged tracks having an 
impact parameter in the plane transverse to the beam less than 0.5 cm from the beam-line. The 
event must have a total measured energy in the laboratory greater than 4.5 GeV within the fiducial 
regions for charged tracks and neutral clusters with energy above 30 MeV. To help reject continuum 
background, the ratio of the second-to-zeroth Fox- Wolfram moment [1C] must be less than 0.5. 

Only upper limits on these modes have been reported and therefore the selection criteria are 
optimized, prior to looking at the signal yield in data, to maximize es/\/ n B- The signal efficiency 
es is evaluated from Monte Carlo, while the background ns is evaluated with data sidebands and 
confirmed by Monte Carlo simulation. The criteria are then made uniform among decay modes 
when appropriate. 

The decays B° -> D [ s *^ tt are reconstructed in the modes D* + — > D^j, Df — > 0vr + , K®K + 
and K*°K + , with </> -> K+K~, K° s vr+vr- and K*° -> K~tt + . The K + and tt - track 
candidates are required to originate from a vertex consistent with the e + e~ interaction point. 
The K® candidates are reconstructed from two oppositely-charged tracks with an invariant mass 
493 < M n + n - < 501 MeV/ c 2 . The 4> candidates are reconstructed from two oppositely-charged 
kaons with an invariant mass 1009 < M K+K - < 1029 MeV/ c 2 . A track is identified as a kaon us- 
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ing the information from the energy-loss measurement in the SVT (only for tracks for momentum 
p < 0.5GeV/c), the DCH (if p < 0.6GeV/c) and the measured Cherenkov angle in the DIRC (if 
p > 0.6GeV/c). Two likelihood selections based on these quantities are used in this analysis: either 
very loose criteria with 95% efficiency and 20% pion contamination, or tight criteria with 85% effi- 
ciency and 5% pion contamination. If not otherwise specified the former is adopted. To reconstruct 
K*° candidates, a K + and a tt~ are required to have an invariant mass 856 < M K + 7T - < 936MeV/c 2 . 
The polarization of the K*° and <p mesons in the Df decays is also utilized to reject backgrounds 
through the use of the helicity angle Oh, defined as the angle between one of the decay products 
and the direction of flight of the meson, in the meson rest frame. Background events are distributed 
uniformly in | cos Oh \ , while signal events are distributed as | cos H \ 2 - K*° candidates are therefore 
required to have | cos0h\ > 0.4, while for the <f> candidates we require | cos6*#| > 0.5. In order to 



reject background from D + 



K° s K+ or K* 



or K*° 7r + , the additional kaon in the reconstruction of Df 
K + is required to pass the tight criteria. Finally, the Df candidates are required 
to have an invariant mass within 10 MeV/c 2 of the nominal mass 

D* + candidates are reconstructed by combining and photon candidates. The selection 
is re-optimized for the B° — > D* + ir~ case, but the resulting selection criteria are very similar to 
that of the B° — > Df-ir~ case. Photons that form a tt° candidate in combination with any other 
photon with energy greater than 70 MeV are not considered. The mass difference between the D* + 



and the the candidate is required to be within 14 MeV/c of the nominal value [ 11 1 
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Figure 1: Distribution of the Fisher discriminant T in B — > Dfir~ Monte Carlo, B — > D*~n + 
data, which has the same form as the signal, and AE* sideband data, which has the same form as 
the background. 



B meson candidates are reconstructed from the Df or D* + candidates and a charged track, 
which is required to fail the tight kaon criteria. Since the B mesons are produced via e + e~ — > 
T(4S*) — > BB, the energy of the B meson in the center-of-mass frame is the beam energy, E£ . 
The distribution of AE* = E* B — Et peaks at zero for the signal. The AE* signal band is 
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defined by |Ai£*| < 36MeV while the AE* sidebands are defined as the remaining events with 
\AE*\ < 120 MeV. 

Backgrounds coming from other B decays, such as B — > D~ir + ,p + with D~ — ► Kgir~ or 
K*°ir~ deserve particular attention because their reconstructed B mass has a shape very similar 
to that of the signal (peaking background). In these events one of the pions in the D~ decay is 
misidentified as a kaon and the distribution of the reconstructed invariant mass of the D~ decay 
products overlaps with the signal Df mass distribution. Nonetheless, since their AE* distribution 
does not peak at AE* = 0, their contribution in the signal region is small. Modes with the same 
final state as the signal but with no intermediate Df can also constitute a peaking background. 
They have the same AE* distribution as the signal, but not the same distribution of the invariant 
mass for the Df candidate. 

In order to reject events where the Df comes from a B candidate and the pion from the other 
B, we require the two candidates to have a probability greater than 0.25% of originating from the 
same vertex, which is 98% efficient on the signal. The remaining background is predominantly 
from continuum qq production with a Df, a or a K*° meson produced in the hadronization of 
one of the two quarks. We use event topology differences between signal and background to reduce 
the continuum contribution. We compute one thrust axis using only the B meson decay product 
candidates and one with all the other tracks. The angle between the two thrust axes (6t) is used to 
discriminate the background. In the center-of-mass frame, BB pairs are produced approximately 
at rest and produce a uniform | cos 9t\ distribution. In contrast, qq pairs are produced back-to-back 
in the center-of-mass frame, which results in a |cos#t| distribution peaking at 1. Depending on 
the background level of each mode we require either | cos#r| < 0.8 or < 0.7. We further suppress 
backgrounds using a Fisher discriminant T constructed from the scalar sum of the center-of-mass 
momenta of all tracks and photons (excluding the B candidate decay products) flowing into 9 
concentric cones centered on the thrust axis of the B candidate (T^j. The more spherical the event, 
the lower the value of T . Figure |l| shows the distribution of this variable in data sidebands, which 
have the characteristics of the background, in simulated signal events and in a control sample of 
approximately 1500 B° — > D*~ir + fully reconstructed events with D*~ — > D°ir~ and D° — ► K + ir~ . 
It is a copius decay channel, with low background, similar final state, therefore well suited to 
investigate signal properties in data. The signal distribution is well reproduced by the simulation. 
We require T < 0.05 for the Df -» 0vr+ and K° K+ modes and T < 0.2 for Df -» K*° K+. 
The selection criteria for B° — ► D* + 7r~ are optimized separately, but are close to the ones for 
B° -> Df-R-. 

As a measure of the B meson mass, the beam-energy substituted mass is defined as toes = 



Efe am — p*g, where pg is the momentum vector of the B meson candidate in the center-of- 
mass frame, calculated from the measured momenta of the decay products. The toes distribu- 
tion for the signal is well described by a Gaussian distribution dominated by the resolution of 
the beam energy measurement, and therefore independent of the decay mode. The combinato- 
rial background is empirically described by a threshold function (the so called "ARGUS" shape 



dN/dniES = ABrriEsyl - mES 2 /E^ am exp [-£ (l - mES 2 /E^ am )] , a function introduced by the 
ARGUS Collaboration |l3j ) for each mode. 

Figure [2] shows the toes distribution for each of the modes. An unbinned maximum-likelihood 
fit is used to fit the toes distributions for signal and ARGUS shaped background contributions. The 
different Df decay modes are combined. The mean and width of the signal distribution are fixed to 
the values obtained in a high statistics B° — ► D^*^tt + sample. The ARGUS shape parameter [13| 
is fixed to the value fitted on the data after having released, in order to increase the statistics, the 
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Figure 2: Distribution of mES for the B° — > D+tt~ (left) and B° — > D* + ir~ (right) candidates after 
all selection criteria, including the AE* signal window. The fits used to obtain the signal yield are 
described in the text. The contributions from the individual modes are also shown. 



requirements on the Df mass (accepted within 40MeV/c 2 of the nominal mass) and the A.E*(to 
50MeV). The signal yields and combinatorial background, as returned from the fit, are given in 
Table § 

Table 1: The number of events in the signal box (N s i g i, ox ), the signal yield (N s i g ) and the com- 
binatorial background (N com b) as extracted from the likelihood fit, the efficiency (e), the peaking 
background (Np^), and the measured branching fraction (£>). N s i g and ^V com 6 and B are not 
available for modes with too few events. N pea k is not reported if no event is found in the Df mass 
sideband. Only statistical errors are quoted. 



Df Mode 


£ 


Nsigbox 


N sig 


Ncomb 


Npeak 


B 




% 


events 


events 


events 


events 


10- 5 






16.6 


7 


6.4 ± 2.7 


1.2±0.5 




3.5 ±1.5 


K*°K+ 


9.7 


6 


5.1±2.4 


1.9±0.6 


2.3 ± 1.8 


2.2 ± 1.8 


K° S K+ 


12.2 


4 


3.4±2.0 


1.2±0.5 




3.7 ± 2.2 


total 




17 


14.9±4.1 


4.4±0.9 


2.3 ± 1.8 


3.1 ±1.0 


B° -> D-+-K- 




7.8 


1 










K*°K+ 


3.3 


3 


2.9 ± 1.8 


0.17±0.17 


0.20 ± 0.14 


6 5+ 36 

D.O_ 2 g 


K°K+ 


5.1 













total 




4 


3.5 ± 2.0 


0.94±0.38 


0.20 ± 0.14 


o i+U.8 
£ - L -l.Q 



The efficiency for the selection requirements is given in Table |j] as obtained from simulation. 
The measurement of the branching fractions for the individual modes are shown as a cross check 
but they are not used to obtain the result. The branching fraction is determined from the signal 
yield (N s i g ), the peaking background (N pea k), the efficiency (e) and the total number of BB events 
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in the sample (Nbb = 61.6 ± 0.6 x 10 6 ). The fit separates the combinatoric background (N com i,) 
from N s i g which is the sum of the signal and the peaking background. We estimate the total 
peaking background by fitting the tties distribution in the Df mass sidebands. The observed yield 
is rescaled to the signal region based on the Df mass distribution in background events. The 
only peaking background that would not be included in this estimate comes from B decays with 
a in the final state. Simulation of a large number of events in these decay modes shows their 
contribution to be negligible. 

Within a 3a uies window we find 17 B° — > Dfir~ and 4 B° — ► D* + -k~ events, and a gaussian 
component (signal and peaking background) of 14.9±4.1 and 3.5±2.0 events, respectively. The 
B° — ► Dfir~ yield has a probability of 4.4 x 10~ 4 to be a fluctuation of the background (3.5a) 
and we measure a branching fraction B(B° — ► Dfn~) = (3.1 ± 1.0(stat.)) xl0~ 5 . The B° — > 
.D* + 7r~ yield has a probability of 2.2<t to be a background fluctuation and the measured branching 
fraction is B(B° — » D* + 7r~) = (2.ll5o (stat.)) xlCT 5 . Probabilities are computed including all the 
uncertainties on the backgrounds. 

As a consistency check for the B° — ► Dfir~ search, we plot in Figure |3] the AE* projection 
for the D+tt~ and the D* + tt~ modes after requiring mES he within a 3a window of the known B 
mass. A comparison of the observed AE* distribution with the expectations from the combinatorial 
background, the component of the peaking background from B° — ► D~n + and B° — > D~p + , and 
the signal itself shows good agreement. The other peaking components with the same final state 
as the signal are included in the signal. 




-0.2 -0.1 0.1 0.2 -0.2 -0.1 0.1 0.2 

AE*(GeV) AE*(GeV) 



Figure 3: AE* projection for B° — ► Dfir~ (left) and B — > .D* + 7r~ (right) events in a 3<7 tties 
window around the B mass. The signal from the Monte Carlo simulation (open histogram), the 
combinatorial background (cross hatched histogram), and the B° — > Z)~7r + and 5° — ► D~p + 
components (hatched histogram) are overlaid. 

The total systematic error is the sum in quadrature of the contributions shown in Table |2|. 
The systematic uncertainty on the combinatorial background subtraction derives from varying the 
background ARGUS shape within the statistical uncertainty on its determination. The uncertainty 
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on the peaking background accounts for the limited size of the sample used to estimate it. The 
uncertainties on the Df and D* + branching fractions are taken from [11]: the dominant source 
is B(Dg — > <J)tt + ), a 25% relative error correlated among all modes, since the other branching 
fractions are measured relative to it. The uncertainty due to the possibility that the simulation does 
not appropriately reproduce the shape of the event selection variables is estimated by comparing 
the corresponding distributions between signal simulation and a copious and pure B° -> D*~ 
control sample. The tracking efficiency is computed from a sample of e + e~ — ► t + t~ events, 
with one r decaying into three tracks and one neutrino, and one decaying into one track and 
one neutrino. We estimate the K® efficiency uncertainty by comparing the momentum and flight- 
distance distributions in data and Monte Carlo simulation. The kaon identification efficiency is 
derived from a sample of D* + — ► D°tt + ,D° — ► K~ir + decays. 



Table 2: The systematic uncertainties in the measurement of the branching fraction B. 



Uncertainty (xlO 5 ) 





D+7T- 


D*+ vr- 


and D* + branching fractions 


0.82 


0.53 


Peaking Background 


0.44 


0.09 


Selection variables 


0.30 


0.18 


Tracking and K® efficiency 


0.17 


0.11 


Kaon identification 


0.14 


0.09 


Combinatoric background 


0.09 


0.06 


Simulation statistics 


0.05 


0.07 


N B b 


0.04 


0.02 


Photon efficiency 




0.03 


Total 


1.01 


0.59 



In conclusion, we observe 17 B° — ► Dfir~ and 4 B° — ► D* + ir~ candidates in the signal re- 
gion. We therefore report a 3.5<r signal for the b —> u transition B° — > Dfir~, with B(B° — ► 
Dfir~) = (3.1 ± 1.0 (stat.)rfcl.O (syst.)) xlO -5 . Given that the dominant uncertainty comes from 
the knowledge of the Df branching fractions we also compute B(B° — ► Dfir~) x B(Df — > (j)ir + ) = 
(1.11 ± 0.37 ± 0.22) x 10" 6 . The search for B° D* + Tr~ yields a result of 2.2cj significance, and 
a 68% confidence interval B(B° -> D* + it-) = (2.1^* (stat.)±0.6 (syst.)) xl0~ 5 , which can be 
translated into an upper limit B(B° -> D* + tt") < 4.3 x 10~ 5 at 90% C.L. 
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